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Abstract 

Flows through three reference nozzles have been calcu- 
lated to determine the capabilities and limitations of the 
widely used Navier-Stokes solver, PARC. The nozzles 
examined have similar dominant flow characteristics as 
those considered for supersonic transport programs. Flows 
from an inverted velocity profile (IVP) nozzle, an 
underexpanded nozzle, and an ejector nozzle were exam- 
ined. PARC calculations were obtained with its standard 
algebraic turbulence model, Thomas, and the two-equation 
turbulence model, Chien k-e. The Thomas model was run 
with the default coefficient of mixing set both at 0.09 and a 
larger value of 0.13 to improve the mixing prediction. 
Calculations using the default value substantially 
underpredicted the mixing for all three flows. The calcula- 
tions obtained with the higher mixing coefficient better 
predicted mixing in the IVP and underexpanded nozzle 
flows but adversely affected PARC’s convergence charac- 
teristics for the IVP nozzle case. The ejector nozzle case did 
not converge with the Thomas model and the higher mixing 
coefficient. The Chien k-e results were in better agreement 
with the experimental data overall than were those of the 
Thomas run with the default mixing coefficient, but the 
default boundary conditions for k and £ underestimated the 
levels of mixing near the nozzle exits. 

Nomenclature 

A cross-sectional area in ejector nozzle mixing region 
A + van Driest damping constant = 26 

C el Chien k-e turbulence model constant = 1.35 

C e2 Chien k-e turbulence model constant =1.8 

C^ Chien k-e turbulence model constant = 0.09 

f^, f p f 2 terms in Chien k-e turbulence model 
H distance from centerline to top or bottom wall of 

ejector nozzle 

K Von Karman constant = 0.4 1 

k turbulent kinetic energy 

t turbulent mixing length 

£ 0 Thomas model constant (default value = 0.09) 

Re t Reynolds number based on turbulent quantities 
t time 
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u,v velocities 

x,y cartesian coordinates 

y + vertical distance normalized with shear length scale 

P ejector nozzle mixing parameter 

e rate of turbulent kinetic energy dissipation 

|i dynamic viscosity 

|i t turbulent viscosity 

n production term in Chien k-e model 

p density 

o k Chien k— e turbulence model constant =1.0 

G e Chien k-e turbulence model constant =1.3 

co vorticity 

co c Thomas model maximum vorticity 
Subscripts: 

i, j computational coordinates 

max maximum 

min minimum 

Introduction 

In the High-Speed Research (HSR) program, NASA and 
industry are using computational fluid dynamics (CFD) 
codes to analyze flow fields from low-noise nozzle concepts 
that could be installed on a High-Speed Civil Transport 
(HSCT). Supersonic transport development programs of 
previous decades, such as the supersonic transport (SST) 
program, 1 relied almost exclusively upon wind tunnel tests 
to evaluate propulsion system components. Recent ad- 
vances in computer technology and flow-solving methods 
have made CFD codes a useful complement to wind tunnel 
tests. 

CFD has been successfully employed to predict the inter- 
nal flow quantities such as thrust coefficients, discharge 
coefficients, and pressure distributions, but only for nozzles 
with simple geometries, which produce single streams with 
no regions of separations. 2 " 4 Nozzle flows involving the 
mixing of a primary stream with lower energy secondary 
streams are much more difficult to calculate. These flows 
are, however, characteristic of those considered for HSCT 
application. The mixer-ejector nozzles of the current HSR 
program, for example, are designed to entrain ambient air 
that mixes with the high-energy core flow. The reduced 
maximum velocities and temperatures of the mixed flow at 
the nozzle exit plane produce less noise than an unmixed 


nozzle. Recent applications of Full Navier-Stokes (FNS) 
analysis to such mixer-ejector nozzles are presented in Refs. 

5 and 6. 

Because FNS codes are now being relied upon to analyze 
complex HSR nozzles, it is important to determine the 
capabilities of currently available codes and to identify 
aspects that require improvement. As a step in this direction, 
the multipurpose FNS solver PARC is applied to the follow- 
ing three reference nozzle flows: (1) an inverted velocity 
profile nozzle, (2) an underexpanded nozzle with the jet 
exiting into Mach 0.6 external flow, and (3) an ejector 
nozzle. These three cases each involve the mixing of a high- 
energy nozzle stream with a secondary stream (or streams) 
of lower energy, and are characteristic of nozzles currently 
considered for use in the HSR program or previously consid- 
ered in past supersonic transport programs. They provide 
data sets appropriate for evaluating the capability of codes 
like PARC to predict quantities indicative of mixing (veloci- 
ties and temperatures). These cases, however, require sig- 
nificantly less grid generation effort and CPU time to obtain 
a converged solution than the actual nozzles being consid- 
ered in the HSR program (such as those described in Refs. 5 
and 6). 

The following discussion begins with a brief description 
of the code used in this study, PARC. Then, for each of the 
three test cases, the experimental background is presented 
followed by a description of the procedure employed to 
obtain calculations and the comparison with experimental 
data. Effects of turbulence model selection on flow field 
predictions are emphasized in the discussion. 

The PARC Code 

The PARC code 7 ’ 8 is an internal flow Navier-Stokes code 
used extensively by government and industry to analyze 
propulsion flows. PARC was derived from the ARC exter- 
nal flow Navier-Stokes code. 9 ’ 10 One version of the PARC 
code contains the two-dimensional and axisymmetric solver 
(PARC2D), whereas the other version contains the three- 
dimensional solver (PARC3D). The governing equations of 
motion are the time-dependent, Reynolds-averaged Navier- 
Stokes equations satisfying a perfect gas relationship and 
Fourier’s heat conduction law. These equations are 
expressed in conservation law form with respect to general 
curvilinear coordinates and are solved with the Beam and 
Warming approximate factorization algorithm. 1 1 A time- 
dependent solver based on a multistage algorithm of Jameson 
has been incorporated into PARC. 12 PARC, however, is 
generally used for steady-state flow simulations. 

Both algebraic and two-equation turbulence models are 
currently available in PARC to analyze turbulent flows. The 
standard algebraic turbulence model is based on the work of 
P.D. Thomas. 13 This model calculates turbulent viscosity 


near surfaces (wall-bounded part of the model) and in 
regions where two or more flows are mixing (free-shear 
layer part of the model) but was optimized for the latter. The 
Baldwin-Lomax algebraic turbulence model 14 is also avail- 
able in PARC. Baldwin-Lomax only calculates turbulent 
viscosity in wall-bounded regions. 

Algebraic turbulence models often model complex flow 
cases inadequately because the single mixing length distri- 
butions used to calculate turbulent viscosity are normally 
tuned to a particular case and are not applicable to all flows. 
Two-equation models avoid this single mixing length limi- 
tation by solving two additional transport equations to calcu- 
late turbulent viscosity, but they are more computationally 
expensive. The Chien low Reynolds number k-e model 15 
with modifications for compressibility added by Nichols 16 
is available in the two-dimensional, axisymmetric code 
(PARC2D). 

Two turbulence models were used with PARC2D in this 
nozzle study to obtain calculations for the three flow cases: 
the Thomas model (the default algebraic turbulence model) 
and the Chien k-e model. The details of these two models 
are presented next. 


Turbulence Models in PARC 
1 . Thomas Model 

The Thomas algebraic turbulence model installed in PARC 
uses Prandtl’s mixing length hypothesis to calculate turbu- 
lent viscosity 

Ht = P^ 2 M 0) 

The turbulent viscosities for wall-bounded regions and 
unbounded regions are calculated separately, based on the 
length scale distributions. For wall-bounded regions the 
turbulent mixing length is described as 

< = Ky(l - e" y+/A+ ) (2) 

For unbounded regions the turbulent mixing length is 

e _ lo [max flujl) - min fluj)] (3) 

COc 

where £ 0 is an adjustable constant that is input as the variable 
COFMIX in the PARC code. The default value of COFMIX 
is 0.09. Because previous flow simulations using the 
Thomas model with i Q = 0.09 have underpredicted the extent 
of mixing, 17 the three nozzle cases of this study were 
examined by using the Thomas model with i Q = 0.13 in 
addition to the default value. This larger value of i Q was 
selected because it provides twice the turbulent viscosity for 
a given density and vorticity field. 
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2. Chien k-£ Model 

The low Reynolds number k-8 model of Chien solves for 
the turbulent viscosity (|i t ) as a function of turbulent kinetic 
energy (k) and rate of turbulent kinetic energy dissipation (£) 
as follows: 

)»it = Cji fjj. p k / £ (4) 


The quantities k and e are solved from the two transport 
equations 
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Discussion of Flow Cases 

Three PARC calculations were run for each of the test 
cases described in following sections. The first used the 
Thomas model with the default i Q = 0.09, and the second 
used Thomas with i Q = 0.13. These will be referred to as 
Thomas (0.09) and Thomas (0. 13) in the rest of this discus- 
sion. The third calculation used the Chien k— £ model. 

1. Inverted Velocity Profile Nozzle 


nozzles maintain an outer nozzle flow at a higher velocity 
than the inner nozzle flow. They have been shown to reduce 
noise relative to a reference conical nozzle. 20,21 

In the current study, an I VP nozzle tested by von Glahn 22,23 
was investigated with PARC. Figure 1 is a schematic of the 
geometry, and Fig. 2 shows the computational grid near the 
nozzle exit plane. The grid was packed to the walls in the 
vertical direction to resolve the developing boundary layers, 
and it was packed to the exit plane in the horizontal direction 
to resolve the sharp gradients of the interacting flows. The 
dark portions of Fig. 2 indicate the locations of high grid 
density. Because the flow was axisymmetric, the bottom 
grid line is the axis of symmetry. The grid had 25 1 horizontal 
points and 191 vertical points, and it extended 1000 cm 
downstream of the nozzle exit plane (approximately 100 
inner nozzle diameters) and 90 cm above the axis of symmetry. 

The flow case investigated with PARC had a primary 
outer flow operating at a total temperature of 1125 K and 
Mach 0.67 and a secondary inner flow operating at a total 
temperature of 290 K and Mach 1 .00, as measured at the 
nozzle exit plane. The outer flow exited the nozzle at a 
higher velocity, but lower Mach number, than the inner flow 
because of its much higher temperature. The surrounding air 
was at rest. The inflow conditions for PARC were con- 
structed by using one-dimensional isentropic relations to 
determine the total conditions that would provide the veloci- 
ties and temperatures measured in the experiment near the 
exit plane. This was done because upstream experimental 
conditions were not specified. 22,23 Typically, nozzle total 
pressures and temperatures are specified at the inflow for 
PARC calculations. Because PARC has difficulty in con- 
verging flows at low speeds (or still air, as in this case), the 
free stream was modeled as having a forward velocity. 

Several difficulties were encountered in the initial 
attempts to obtain solutions for this IVP nozzle case. First, 
the packing of the initial grid was too sparse at the nozzle exit 
plane. Increasing the grid density to that shown in Fig. 2 
helped resolve the sharp flow gradients between the nozzle 
flows and the ambient air. Second, the external flow sepa- 
rated from the aft-facing surface. Moving the free-stream 
inflow boundary to the nozzle exit plane eliminated the area 
of separation and assisted convergence. Finally, because the 
surrounding air in the experiment was at rest, it was desirable 
to run the PARC calculations with as low a free-stream Mach 
number as possible. The initial calculations were run with a 
free-stream Mach number of 0.60. This was gradually 
reduced to the lowest Mach number that provided converged 
solutions for all three PARC calculations, 0.45. 


One of the major environmental challenges NASA and 
industry are addressing in the HSR program is controlling 
the jet noise that will be produced by the HSCTs high thrust 
engines. 18,19 One concept that has been examined exten- 
sively is the inverted velocity profile (IVP) nozzle. IVP 


Velocity profiles obtained from the PARC calculations 
are compared to experimental data in Fig. 3. At a location 
20 cm downstream of the nozzle exit (Fig. 3a), the maximum 
velocities of the inner flow (from the axis of symmetry to a 
radial position of approximately 5 cm) and the outer flow 
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(from 5 cm to approximately 8 cm) match the data well. This 
indicates that the nozzle inflow boundary conditions were 
calculated accurately. Further downstream, at 40 cm 
(Fig. 3b), 81 cm (Fig. 3c), and 130 cm (Fig. 3d), none of the 
three PARC solutions show a flow decay as rapid as the 
experimental data. The Thomas (0.09) solution shows the 
least decay at 81 cm and 130 cm downstream of the nozzle 
exit. The static temperature comparisons in Fig. 4 demon- 
strate the same trends. At 1 30 cm downstream (Fig. 4d), the 
Chien k-e solution shows more temperature decay than the 
other flow solutions but still substantially less than the 
experimental data. 

Although the inability of the three PARC calculations to 
produce the same flow decay as the experimental data was 
largely due to the presence of the Mach 0.45 free-stream 
flow in the calculations, the differences among the solutions 
were due to the significantly different turbulent viscosity 
levels that were calculated for each case. Figure 5 shows 
turbulent viscosity contours (p/|i) for the three PARC solu- 
tions. The Thomas (0.09) solution (Fig. 5a) calculated less 
turbulent viscosity than the other two solutions throughout 
the flow field, which resulted in the poorest velocity and 
temperature decay characteristics. The Thomas (0. 1 3) solu- 
tion (Fig. 5b) shows higher levels of turbulent viscosity than 
the Thomas (0.09) solution, but its contours are discontinu- 
ous. The Thomas (0.13) solution was the only one of the 
three that remained somewhat oscillatory and had some 
difficulty in converging. 

Because the Chien k-e model calculates turbulent viscos- 
ity by solving transport equations, boundary conditions for 
k and e must be specified at each of the boundaries, including 
the inflow. The quantities k and e were not measured in the 
experiment, so the default boundary conditions for the 
inflow, extrapolations of the values for k and £ from the 
interior, were used. Equation (4) shows that the calculated 
turbulent viscosity is proportional to k 2 . Near the inflow 
boundary, the values of k were quite small, and as a result, 
the Chien k-£ solution’s turbulent viscosity contours (Fig. 
5c) are quite low just downstream of the nozzle exit (as low 
or lower than the two Thomas solutions) but much higher 
further downstream. 

2. Underexpanded Nozzle 

The second flow case investigated was the underexpanded 
nozzle tested by Heltsley and Crosswy. 24 Unlike the IVP 
nozzle, the underexpanded nozzle was tested in a wind 
tunnel operating at a Mach number of 0.6; therefore, the free- 
stream Mach number for the PARC calculations was not 
specified artificially. A schematic of this axisymmetric 
single flow nozzle is shown in Fig. 6. The grid near the 
nozzle exit plane is shown in Fig. 7. The internal surface of 
the nozzle upstream of the throat was modified by removing 
a sharp corner to simplify grid generation. This grid had 25 1 
horizontal points and 141 vertical points, and it extended 


60.96 cm in the horizontal direction from the nozzle exit 
plane and 1 5.24 cm in the vertical direction from the axis of 
symmetry (bottom grid line). The flow case investigated had 
the following operating conditions : a wind tunnel stagnation 
pressure of 136.4 kPa, a free-stream Mach number of 0.60, 
a free-stream stagnation temperature of 355.5 K, a nozzle 
Mach number of 1 .563 at the exit plane, a nozzle stagnation 
temperature of 355.5 K, and a nozzle static pressure ratio 
(static pressure of the flow at the nozzle exit plane divided by 
the free-stream static pressure) of 4.5. The three PARC 
calculations obtained for this nozzle flow converged rela- 
tively easily and did not demonstrate the grid or free-stream 
difficulties that the IVP nozzle presented. 

The experimental data indicates that a normal shock 
occurs near x = 2.50 cm. All of the PARC calculations 
predicted the shock to occur at approximately 2.30 cm 
downstream of the nozzle exit, and thus they are in good 
agreement with the experimental data. Adamson and 
Nicholls 25 presented data that indicate an underexpanded 
nozzle operating under similar conditions but exiting into 
still air has a normal shock occurring at approximately 
2.00 cm downstream of the nozzle exit. 

Figure 8 shows velocity profiles at four stations down- 
stream of the nozzle exit. At 1.27 cm downstream of the 
nozzle exit (Fig. 8a), the velocities are much higher than the 
experimental data. The large difference at this location may 
be due in part to aluminum oxide particles, as measured by 
laser Doppler velocimetry (LDV), lagging in velocity rela- 
tive to the highly accelerating flow leaving the nozzle. 
Further downstream, the solutions vary widely. Good agree- 
ment with experimental data is achieved only with the 
Thomas (0.013) solution at 6.35 cm (Fig. 8c) and 8.90 cm 
(Fig. 8d) downstream. The Thomas (0.09) and Chien k-£ 
solutions demonstrate significantly less core flow decay. 

The differences in these velocity profiles relative to the 
experimental data are due to the variation in turbulent 
viscosities calculated in the flows. Turbulent viscosity 
contours are shown in Fig. 9. The Thomas (0.13) solution 
(Fig. 9b) has turbulent viscosity levels that are much higher 
just downstream of the nozzle exit. The Chien k— £ turbulent 
viscosity levels (Fig. 9c) do not become substantially higher 
than the Thomas (0.09) levels (Fig. 9a) until well down- 
stream of the last position where experimental velocities 
were available for comparison (x = 8.90 cm). As for the IVP 
nozzle case, the Chien k-£ calculation used the default 
inflow boundary conditions with the turbulent quantities 
extrapolated from the interior. Using more accurate values 
for k and £ at the inflow may have improved the agreement 
with experimental data. 

3. Ejector Nozzle 

The last test case examined was a two-dimensional ejector 
nozzle, which is a simplification of the mixer-ejector nozzles 


4 


considered in the HSR program. Mixer-ejector nozzles 26 
entrain large amounts of secondary air that mix with the high 
energy flow from the engine core to lower the average jet 
velocity and resultant jet noise. High thrust levels are 
maintained despite the lower average jet velocity due to the 
mass augmentation of the secondary airflow. Mixer-ejector 
nozzles are a specific class of ejector nozzles that utilize 
devices such as lobed chutes to enhance mixing. 

Ejector nozzles have received much attention for their 
potential application in providing thrust augmentation for 
vertical and short takeoff or landing (VSTOL) aircraft. 
Bevilaqua discusses some of the fundamentals of ejector 
nozzle operation and performance in Refs. 27 and 28. 
Analyses of ejector nozzle concepts began with, and contin- 
ues to include, theoretical treatments. In recent years, CFD 
has also been applied to ejector nozzle concepts. 5,6 ’ 29 

For the current study, a two-dimensional ejector nozzle 
tested by Gilbert and Hill 30 was investigated with PARC. 
Figure 10 shows the geometry of the ejector test rig used in 
their experiments. The flow case that was investigated with 
PARC is referred to as run number 9 in Ref. 30. This flow 
case had the following operating conditions: a primary 
nozzle total pressure of 246 kPa, a primary total temperature 
of 358 K, an ambient pressure of 101 kPa, and an ambient 
temperature of 305 K. Figure 1 1 shows the details of the 
computational grid used in the PARC calculations in the 
vicinity of the primary nozzle and secondary flow inlet. The 
upstream portion of the secondary flow inlet was modeled as 
shown in Fig. 1 0 to avoid grid skewness problems that would 
have occurred with exactly modeling the ejector wall far 
upstream. The grid had 25 1 horizontal and 121 vertical grid 
points, and it extended to 26.7 cm downstream of the nozzle 
exit plane (immediately before the diffuser entrance). The 
x-axis is a line of symmetry, and positions in the mixing 
section are measured relative to the primary nozzle exit 
plane (x = 0 cm). 

In obtaining calculations with PARC, only the Thomas 
(0.09) and Chien k-e models provided converged solutions. 
The Thomas (0.13) calculation’s large, but discontinuous, 
turbulent viscosities caused sporadic mixing through the 
ejector nozzle that was not compatible with the free bound- 
ary condition at the outflow (pressure and temperature are 
specified). The Thomas (0. 13) model calculated turbulent 
viscosities that were discontinuous for the IVP and under- 
expanded nozzles, but for those cases, the free boundary 
condition was specified far downstream and did not conflict 
with the mixing flows. As a result, only the Thomas (0.09) 
and Chien k-e solutions are compared for this ejector nozzle 
case. 

Figure 12 shows velocity profiles at four stations in the 
mixing section. The PARC calculations were obtained for 
only half of the test section because of symmetry. The 
experimental data, obtained for the entire height of the test 


section, showed some asymmetry which may be observed in 
Fig. 12, particularly near the mixing region centerline. For 
each axial station, the velocity profiles are plotted versus 
normalized vertical position (y/H), where H is the local 
distance from the centerline to either the top wall or bottom 
wall. These velocity profiles show that the Chien k-e 
solution matches the experimental velocity profiles much 
more closely than the Thomas (0.09) solution. The Chien 
k-e solution shows that the mixing close to the nozzle exit 
is underpredicted somewhat, which again may be due to the 
inflow boundary conditions used for k-e. The agreement 
with experimental data improves further downstream of the 
nozzle exit. 

A measure of mixing effectiveness proposed by 
Bevilaqua 28 relates the mixing effectiveness ((3) to the flat- 
ness of the velocity profile: 

f v 2 dA 

( 12 ) 

v 2 A 

This parameter was calculated for the PARC solutions and 
the experimental data. Figure 13 compares the mixing 
effectiveness obtained from the PARC calculations using 
both turbulence models with the experimental data and 
shows that the Chien k-e solution agrees best with the 
experimental data. Figure 14 shows total temperature pro- 
files at two positions in the mixing section. Neither solution 
predicts the correct decay of the temperature in the center of 
the flow. Throughout the mixing section, the k-e solution 
agrees much better with the experimental data than does the 
Thomas (0.09) solution. 

Figure 15 shows turbulent viscosity contours for the 
Thomas (0.09) and Chien k-e solutions. The Chien k-e 
model’s much higher levels of turbulent viscosity are 
directly responsible for better mixing predictions relative to 
the Thomas (0.09) calculation. Although the Chien k-e 
model predicted the mixing behavior quite well, the agree- 
ment with experimental data possibly could have been 
improved if more appropriate boundary conditions for k and 
e were used at the inflow. 

A final comparison of the solutions to experimental data 
is provided in terms of predicted mass flow. Since mass 
augmenting ejectors have been designed to entrain as much 
secondary flow as possible, the pumping ratio (defined as the 
ratio of entrained secondary mass flow to primary mass 
flow) has been used as a measure of ejector effectiveness. In 
the Gilbert and Hill experiment, the pumping ratio was 4.67, 
whereas the PARC solution using the Chien k-e model 
provided a pumping ratio of 4.28 and that using Thomas 
(0.09) predicted 4.17. The CFD results indicate that the 
secondary flow entrainment is relatively independent of the 
extent of mixing. 
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Conclusions 

Navier-Stokes calculations have been compared to 
experimental data for three reference nozzles that have many 
of the same dominant flow characteristics as those consid- 
ered for supersonic transport application. The calculations 
were obtained with the general purpose Navier-Stokes solver 
PARC, using the code’s standard Thomas algebraic turbu- 
lence model (with the default coefficient of mixing, £ Q = 

0.09, and a larger value, £ Q = 0. 1 3, in an attempt to improve 
mixing in the calculations) and the Chien two-equation k— 8 
model. 

The solutions obtained with Thomas (0.09) substantially 
underpredicted the extent of mixing for all three flows. 
Those obtained with Thomas (0.13) provided better results 
for the IVP nozzle and underexpanded nozzle flows relative 
to those obtained with Thomas (0.09), but the larger magni- 
tudes of turbulent viscosity also caused the calculations to 
have poorer convergence characteristics for the IVP nozzle 
case. For the ejector nozzle flow, the Thomas (0. 1 3) calcu- 
lation did not converge. Well downstream of the nozzle exit, 
the Chien k-e model produced somewhat better mixing 
characteristics for the IVP nozzle than did the Thomas 
(0.09), and it agreed much more closely with experimental 
data for the ejector nozzle case. The Chien k-e solution did 
not provide better agreement with experimental data than 
the Thomas solutions did for the underexpanded nozzle 
flow. However, the velocities obtained with LDV were only 
available for axial positions very close to the nozzle exit, and 
it was determined that the inflow boundary conditions used 
for k and e resulted in low levels of turbulent mixing close 
to the nozzle exits of the three cases examined. 

The results of this study indicate that the selection of 
turbulence model has a significant effect on the calculation 
of high-speed nozzle flows. The Thomas model using £ Q = 
0.09 (which was intended for free shear layer calculations) 
appears to be inadequate for the calculation of nozzle flows 
such as those considered here. Using £ Q = 0. 1 3 improved the 
level of mixing for the unbounded free shear layer flows but 
adversely affected PARC’ s convergence characteristics. The 
Chien k-£ model produces better flow field predictions 
overall than Thomas (0.09) but its default inflow boundary 
conditions for k and £ caused an underestimation of mixing 
near the nozzle exit. Using more appropriate boundary 
conditions may help solve this problem. 

Improving the capabilities of Navier-Stokes codes such as 
PARC to calculate complex nozzle flows will require con- 
tinued development and implementation of appropriate tur- 
bulence models (and boundary conditions for turbulence 
models that solve transport equations like k-£). In addition, 
the ability of compressible flow codes to handle static or 


very slow flow conditions (like the air surrounding the IVP 
nozzle) must be improved. The nozzle test cases examined 
in this study can be used to show where progress is being 
made in these areas and to what extent Navier-Stokes solvers 
can be expected to provide reliable nozzle flow predictions. 
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Figure 1. Schematic of Inverted Velocity Profile Nozzle. 



Figure 2. Grid detail near inverted velocity profile nozzle exit plane. 
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Figure 3. Velocity profiles for IVP nozzle. 
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Figure 4. Temperature profiles for IVP nozzle. 
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Figure 5. Turbulent viscosity contours for inverted velocity profile nozzle. 
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Figure 6. Schematic of underexpanded nozzle. 



Figure 7. Grid detail near underexpanded nozzle exit plane. 


12 


RADIAL POSITION (cm) RADIAL POSITION (cm) 




(a) x = 1.27 cm 


(b) x = 3.80 cm 




(c) x = 6.35 cm 


(d) x = 8.90 cm 


Figure 8. Velocity profiles for underexpanded nozzle. 
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Figure 9. Turbulent viscosity contours for underexpanded nozzle. 
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Nozzle exit height = 0.309 cm 
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Figure 10. Schematic of 2D ejector nozzle. 



Figure 11. Grid detail near 2D ejector nozzle exit plane. 
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Figure 1 2. Velocity profiles for 2D ejector nozzle. 
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Figure 1 3. Mixing effectiveness for 2D ejector nozzle. 
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Figure 1 4. Total temperature profiles for 2D ejector nozzle. 
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Figure 15. Turbulent viscosity contours for 2D ejector nozzle. 
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